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Abstract 

We present some new problems in spectral optimization. The first one consists in 

determining the best domain for the Dirichlet energy (or for the first eigenvalue) of 

the metric Laplacian, and we consider in particular Riemannian or Finsler manifolds, 

^^ Carnot-Caratheodory spaces, Gaussian spaces. The second one deals with the optimal 

^Sj shape of a graph when the minimization cost is of spectral type. The third one is the 

^ I optimization problem for a Schrodinger potential in suitable classes. 

-^ Keywords: shape optimization, eigenvalues, Sobolev spaces, metric spaces, optimal 

graphs, optimal potentials. 
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O 1 Introduction 

q 

r^ Spectral optimization theory goes back to 1877, when Lord Raileigh conjectured, in his book 

"^ "The Theory of Sound" [12], that among all drums of prescribed area the circular one had 

ri the lowest sound. Here are his precise words: 

I— —I "// the area of a membrane he given, there must evidently he some form of boundary for 

^^ which the pitch (of the principal tone) is the gravest possible, and this form can be no other 

j> than the circle. . . " 

On Since then, many other optimization problems involving the spectrum of the Laplace 

^ operator have been considered (see for instance the survey paper [5] and the books [3], 

■^ [H] , [E] ) , showing the existence of optimal shapes and their qualitative properties together 

.^ with the corresponding necessary conditions of optimality. However, in spite of the strong 

C^ development of the theory, many problems still remain open and many conjectures are still 

^_j waiting for a proof. 

• • In this paper we present some different directions of research; our goal is to consider 

. ^H spectral optimization issues for the following three classes of problems. 

X 

H • Optimization with respect to the domain for functionals like the Dirichlet energy or the 

first Dirichlet eigenvalue related to the metric Laplacian. This operator is in general 
non-linear and acts on functions defined on a general metric space; of particular interest 
are the cases when the metric space consists in a Riemannian or Finsler manifold, in 
a Carnot-Caratheodory space, in a Gaussian space. 

• Optimization of the shape of a graph with respect to the Dirichlet energy or to the 
first eigenvalue. In this case some explicit examples can be provided, together with 
some general necessary conditions of optimality. 



• 



Optimization of the potential V{x) in a Schrodinger equation of the form — Aii + 
V{x)u = f{x). The potential will be submitted to some suitable integral constraints 
and an existence result will be provided for several cost functionals. 



The three cases above wih be treated in Sections [2| [3] and |4j respectively. In ah the 
cases Dirichlet boundary conditions will be considered; other kinds of boundary conditions 
would require completely different mathematical tools that in many cases are only partially 
developed. Our main concern is addressed to the existence of optimal solutions; other 
very interesting questions, like for instance the regularity of optimal solutions, have at 
present only limited and partial answers. In all the three cases, the existence of an optimal 
domain is obtained through the direct methods of the calculus of variations, that require 
two main ingredients: compactness of the space of competitors and semi-continuity of the 
cost functional. In the literature (see for instance [3j) some useful topologies on the family 
of admissible domains have been introduced, in order to provide the necessary compactness 
properties. The semi-continuity of the cost functional is a more involved issue and requires 
some careful analysis. 

The purpose of the present paper is not to provide new proofs or new results but mainly 
to illustrate the field of spectral optimization problems through some examples and to 
discuss some crucial issues by proposing some interesting problems that, to the best of our 
knowledge, are still open. 

In Section [2] we consider the general framework of metric spaces, on which the metric 
Laplacian operator can be defined, together with the related energy and spectral eigenvalues. 
We recall a general existence result of an optimal domain, obtained in ^ , and we show some 
related examples concerning Riemannian or Finsler manifolds, Carnot-Caratheodory spaces, 
Gaussian spaces. 

In Section |3] we consider the case of spectral optimization problems for graphs, and in 
some cases we are able to provide explicitly the optimal shapes. We consider a natural 
convergence on the set of metric graphs in terms of the connectivity matrices of the graphs 
and the lengths of the edges. It is not hard to check that the spectral functionals we 
consider are continuous with respect to this convergence. On the other hand the family of 
admissible graphs endowed with such a convergence is not even complete, which gives raise 
to some counterexamples to the existence. Thus, we investigate the problem in a wider, 
more appropriate class of competitors. 

In the last Section [4] we consider potentials for Schrodinger equations and the related 
optimization problems. In this case the admissible set of choices is just L^{Q,), the set of 
positive integrable functions on 0, and the constraints are given by some integral inequalities. 
In this case, both the compactness of the optimizing sequences and the semi-continuity of 
the cost functional are quite involved questions, and the existence of optimal potentials is 
only known in some particular cases, leaving several interesting problems still open. 

2 Spectral optimization in metric spaces 

In this section we consider spectral optimization problems in the class of subsets of some 
ambient metric space {X, d) endowed with a finite Borel measure m. We do not assume any 
compactness or boundedness of X with respect to the distance d. Our main assumption is 
the compactness of the inclusion L'^(m) C H^{X,m), where H^{X,m) is a Sobolev space of 
functions on (X,m), which we define in each of the cases we consider. 

2.1 Metric measure spaces 

In [S] we consider a separable metric space {X, d) endowed with a finite Borel measure m 
and a Riesz subspace H of L?'{m) satisfying the Stone property, i.e. 

if u G H, then u Al €z H and \u\ € H. 



Let D : H ^ L'^{m) be a convex, l-homogeneous map which is also local, i.e. 

D{u Vv) = Du- I{u>v} + Dv ■ I{u<v}, Vn, v e H. 

We consider H endowed with the norm 

MH = {\\u\\h + \\Du\\l,f\ 

Moreover, we assume that 

{Til) the inclusion i : H ^^ L^ is compact; 

{712) the norm of the gradient is lower semi-continuous with respect to the L^ convergence, 
i.e. for each sequence Un bounded in H and convergent in the strong L^ norm to a 
function u G L'^{m), we have that u £ H and 



/ \Du\ dm < liminf / |-Dun| dm; 
Jx "^°° Jx 



{T-13) the linear subspace HriC{X), where C{X) denotes the set of real continuous functions 
on X, is dense in H with respect to the norm || • \\h- 

An interesting example of subspace H with the properties above is given by the Sobolev 
space H^ {X, m) in the sense of Cheeger [10] . 
For any set fi C X, we define the space 

Hq{Q) = {u£H : cap({n / 0} \ (7) = O}, 

where the capacity cap(£') of a generic set E C X, is defined by 

cap(£^) = inf {||^t||// : u^H, u>0 on X, u>l in a neighbourhood of £'}. 

Definition 2.1. For each Borel set (7 and each k > 1, we define 

Afc(^) = inf sup { iDu]"^ dm : ue K, u'^dm = l}, (2.1) 

KcHoin) L 7q Jn ^ 

where the infimum is over all k-dimensional linear subspaces K of HqIQ). 

Definition 2.2. For each Borel set Q and each f £ L'^{Q,m), the Dirichlet energy of 0, is 
defined as 

Ef{n)=ini\- / \Du\'^dm + - u^ dm - / uf dm : u£Ho{n)}. (2.2) 

Remark 2.3. In the cases when we have the inequality ||w||j;^2(^) < C||-Dn||j;^2(^), for each 
u £ H, it is more convenient to define the energy Ef{Q) as 

Ef{n) = inf \^ I \Du\'^dm- [ uf dm : ueHo{n)}. (2.3) 

Also in this case the statement of the following theorem remains valid. 

Theorem 2.4. Suppose that {X,d) is a separable metric space with a finite Borel measure 
m and suppose that H C L^{X,m) and D : H ^^ Lp'{X,m) are as above. Then the shape 
optimization problems 

min{Ef{n) : Q C X, m(f]) < l}, 

and 

min{Afe(ri) : QcX, m(0) < l}, 

have solutions, which are quasi-open sets, i.e. level sets of the form {u > 0} for some 
function u £ H. 



Remark 2.5. The existence result of Theorem 2.4 holds, in the same form, for several other 



shape functionals F{Q); the only required assumptions (see j9j) are: 

- F is monotone decreasing with respect to the inclusion, that is 

F(ili) < F(il2) whenever Q2 C ^i; 

- F is 7-lower semi-continuous, that is 

F{Q) < limini F{Qn) whenever wq^ — )> wq in L {X,m) 



where wq is the solution of the minimization problem (2.2) with / = 1 



For instance, the following cases belong to the class above. 

Integral functionals. Given a right-hand side / we consider the PDE formally written as 

—Au + u = f in il., u G Hq{0,), 



whose precise meaning is given through the minimization problem (2.2), and which provides, 
for every admissible domain Q, a unique solution uq that we assume extended by zero outside 
of 0,. The cost F{0,) = J{un) is then obtained by taking 

J{u)= I j(^x,u{x)^dm 
Jx 

for a suitable integrand j. If / > and j(x, •) is decreasing, this cost verifies the conditions 
above. 

Spectral optimization. For every admissible domain we consider the eigenvalues Afc(il) 
and the spectrum X{Q,) = (Afc(il))^. Taking the cost 



of Definition 



2.1 



F{n) = ^{x{n)) 



we have that the assumptions above are satisfied as soon as the function <1> : [0, +00]^ — )• 
[0, -|-oo] is lower semicontinuous and increasing, in the sense that 

A^'^Afc VfcGN ^ $(A) <liminf$(A''), 
Xk<fik Vfc G N => $(A) < $(/x) . 

2.2 Finsler manifolds 

Consider a differentiable manifold M of dimension d endowed with a Finsler structure, i.e. 
with a map F : TM — )• [0, -|-oo) which has the following properties: 

1. F is smooth on TM \ {0}; 

2. F is 1-homogeneous, i.e. F{x,XX) = \X\F{x,X), VA G M; 

3. F is strictly convex, i.e. the Hessian matrix gij{x) = \-qj^jq^[F'^]{x,X) is positive 
definite for each (j;, X) G TM. 



With these properties, the function F{x, •) : TxM — )• [0, +00) is a norm on the tangent 
space TxM, for each x G M. We define the gradient of a function / G C°^(M) as Df{x) := 
F*{x, dfx), where dfx stays for the differential of / at the point x £ M and F*{x, •) : T*M — )• 
M is the co-Finsler metric, defined for every ^ S TjfM as 

F*{x,0= sup -^ 
ydT^M F[x,y) 

The Finsler manifold (M, F) is a metric space with the distance: 

dpix, y) = inf { y F(7(t), 7(t)) dt : 7 : [0, 1] ^ M, 7(0) = x, 7(1) = y}. 

For any finite Borel measure m on M, we define -ff := Hq{M, F, m) as the closure of the set 
of differentiable functions with compact support C'^{M), with respect to the norm 



li/ll-^ -I- II TliiW^ 

I II L^(m) II II L^(m) 



The functionals A^ and Ef are defined as in (2.1) and (2.2), on the class of quasi-open 
sets, related to the H^{M,F,m) capacity. Various choices for the measure m are available, 
according to the nature of the Finsler manifold M. For example, if M is an open subset 
of M , it is natural to consider the Lebesgue measure m = C . In this case, the non- 
linear operator associated to the functional J F*{x, duxY' dx is called Finsler Laplacian. On 
the other hand, for a generic manifold M of dimension d, a canonical choice for m is the 
Busemann-Hausdorff measure mp, i-e. the d-dimensional Hausdorff measure with respect to 
the distance dp- The non-linear operator associated to the functional J F*{x, duxY drnp^x) 
is the generalization of the Laplace-Beltrami operator and its eigenvalues are defined as in 



(2.1). In view of Theorem 2.4, we have the following existence results: 



Theorem 2.6. Given a compact Finsler manifold (M, F) with Busemann-Hausdorjf mea- 
sure mp, the following problems have solutions: 



mm 



mm 



\ Xk{^) '■ mpi^) < c, Q quasi-open, il C M>, 
< Ef{Q) : mp{Q) < c, Q quasi-open, $1 C M >, 



for any k gN, < c< mp{M) and f G L'^{M, mp). 

Theorem 2.7. Consider an open set M C M endowed with a Finsler structure F and the 
Lebesgue measure C^. If the diameter of M with respect to the Finsler metric dp is finite, 
then the following problems have solutions: 



mm 



mm 



< Afc(ri) : \Q\ < c, rj quasi-open, Q, C M>, 
I Ef{Q) : \^\ < c, Q, quasi-open, Q, C M>, 



where k £N, \i}\ denotes the Lebesgue measure ofQ, c is a constant such that < c < \M\ 
and f eL^{M). 

Remark 2.8. In [llj it was shown that if the Finsler metrics F{x, ■) on W^ does not depend 
on X E M , then the solution of the optimization problem 



mm 



|Ai(ri) : \Q\<c,Q quasi-open, f] C M°'|, 



is the ball of measure c. It is clear that it is also the case when in the hypotheses of Theorem 
12.71 one considers c > such that there is a ball of measure c contained in M. On the other 
hand , if c is big enough the solution is not, in general, the geodesic ball in M (see }15j). 
If the Finsler metric is not constant in x, the solution will not be a ball even for small c. 
In this case it is natural to ask whether the optimal set gets close to the geodesic ball as 
c — )• 0. In [18] this problem was discussed in the case when M is a Riemannian manifold. 
The same question for a generic Finsler manifold is still open. 

2.3 Gaussian spaces 

Consider the Euclidean space M^ endowed with the Gaussian measure 

-1 / -^1 "T ■'^2 



m = (27r) exp I I dxidx2. 

Note that an orthonormal basis on L'^{m) is given by the functions Hn^k{xi,X2) := Hn{xi)Hk{x2) 
n,k gN, where Hn : M — t- M are the Hermite polynomials 

(—1)"' 
Hn{x) := ^—I^ exp(xV2) 9^ (exp(-xV2)) , 



which satisfy 

dxHn{x) = ^/nHn-l{x), dlHn{x) - xHn{x) = nHn{x). 

We define the Sobolev space W'^''^{^,m) as 

W^''^{M^,m) = {ueL^{m) : |Vn| E ^^(m)} , (2.4) 

where Vu is the distributional gradient of u. It can be characterized using the basis {//„ ^j^ ^ 
as 

W^'^{R^,m) = {u e L^{m) : ^(n + A;)!*^^^ < +oo}, (2.5) 

n,k 

where Un,k '■= J'^2 Hn^kudm. At this point it is clear that the inclusion W'^''^{^,m) C 
L'^{m) is compact and that the estimate ||m||l2(^) < ||Vn||j;^2(^) holds. Moreover, the linear 
combinations of Hermite polynomials are dense in W^''^{'M?, m) and so C°°(M^)nVF^'^(M^, m) 
is dense in W^''^{M?,m). Thus, we can define the capacity cap(£') of any set E C M? and 
the space Wq' {Q,m) of functions u G VF^'^(M^,m) such that cap({M / 0} n fl'^) = 0. For 



any / S L'^{m), there is a unique w S Wq' {0,,m), which minimizes the functional 



Jf{u) = - / |Vn| dm— / fudm, 
2 Jn Jn 

and defines the energy of il as Ef{i^) := Jf{w). We note that for any v G Wq' (i7,7n) we 
have 



Vw ■ Vv dm = / fw dm, 
n Jn 

and so, we say that w is the weak solution of the problem —Aw + x ■ Vw = / in Wq ' {Q, m) . 
Since ||Vit||x,2(^) < ||/||L2(,m), we have that the operator Rq : L'^{m,) — )• L'^{m), which 
associates to each / G L'^{m) the function Rn{f) '■= w, is compact. Thus Rq is the 
resolvent of an operator — A + x • V, which is the Ornstein-Uhlenbeck operator on il and 
which has a discrete spectrum cr{Q), given by the sequence < Ai(ri) < A2(ri) < . . . . Note 
that, in the case O = M^, the spectrum is given by (t(M^) = {n + k : n. A; G N}. In particular, 
Ai(R2) = and A2(M^) = Xsi^'^) = 1. We also note that the A;-th eigenvalue Xki^) can be 



represented as in (2.1 ) and so, if O 7^ M , then Ai(r2) > 0. Applying Theorem 2.4 we obtain 



the existence of optimal domains for any A^. 



Theorem 2.9. Consider M? endowed with a non-degenerate Gaussian measure m, i.e. with 
invertible covariance matrix. Then, for any k G N, f £ L^(m) and < c < 1, the following 
optimization problems have solutions: 



mm 



mm 



hkin) : JlcM^ m{n) < c\, 



which are quasi-open sets. 



Remark 2.10. Theorem 2.9 also apphes to penahzed problems, i.e. for any A > 0, /c G N 
and f £ L? (m) , there is a solution of the problems 



mm 



mm 



miEf{VL)+Km{Q) : 17 C M^j, 



(2.6) 



(2.7) 



which is a quasi-open set. As we will see in the example below, these problems are sometimes 
easier to threat when comes to regularity questions and qualitative study of the optimal sets. 

Example 2.11. Let / be the constant 1 in M . By Remark 2.10, the problem (2.7) has a 



solution ri, which we assume to be open and with boundary dVL of class C^ (that we expect 
to be true), we can perform the shape derivative of the energy Ei with respect to some 
vector field V regular enough. Indeed, following \\Q\ Chapter 5], let y : M*^ — )■ M'^ be a C^ 
vector field and for each i > small enough, define $t(x) = x -\- tV{x) and ^t = ^t{^)- 
Then, we have 



dEii^t) 



dt 



1 



t=o 



w' dm, 



where w' is the solution of 



-Aw' + X ■ Vw' = 0, in i7. 



1 u;' = —V ■ Vw, on dQ 
We denote with w the (strong) solution of 

—Aw + X ■ Vw = 1, 



(2. 



(2.9) 



w£WQ'^{n,m) 



and integrate by parts in (2.8) obtaining 
dEiiQt) 



dt 



t=Q 



- / {—Aw-\-x-Vw)w'dm = —-— / 
2 Jn 47r Jq^ 



dw 



dn 



V-ne~\''\'^^dn'^~\ (2.10) 



where n is the exterior normal on dO, and w is the energy function on i}, that is the solution 
of the Ornstein-Uhlenbeck PDE 



a, 2/ 



—Aw + X ■ Vw = 1 in O, w £ Wq ' ($7, m) . 

On the other hand, we have 

dm{Qt) 



dt 



= - [ e~\^\"I^V-nd'H^-\ 
=0 27r Jq^ 



(2.11) 



and so, by the optimality of i7, 

(dEi{nt) 



- 



dt 



+ A 



dm{Qt 
dt 



i=0 



for any vector field V. By (2.10) and (2.11) we obtain 

dw 



dn 



'2A 



on dVl. 



Summarizing, we have obtained that if an optimal domain Vt is regular enough, then the 
following overdetermined boundary value problem has a solution: 



—Aw + X ■ 'Vw = 1, in il, 
w = 0, on dQ, 

-\/2A, on dn. 



(2.12) 



dw 
dn 



It is straightforward to check that the following domains satisfy this condition: 

• the half-space Q = {xi > c}, for a given c G M, 

• the strip f] = {|xi| < a}, for some a > 0, 

• the euclidean ball il = {|x| < r}, for some r > 0, 

• the external domain of a ball fi = {|x| > r}, for r > 0. 

We do not know which of these domains is optimal and if there are other domains il for 



which the overdetermined problem (2.12) has a solution 



2.4 Carnot-Caratheodory spaces 

Consider a bounded open and connected set D C M'^ and C°° vector fields Yi, . . . ,Yn defined 
on a neighbourhood U of D. We say that the vector fields satisfy the Hormander's condition 
on U, if the Lie algebra generated by Yi, . . . ,Yn has dimension d in each point x G [/. 

We define the Sobolev space Wq' {D; Y) on D with respect to the family of vector fields 
Y = (Yi, . . . , Yn) as the closure of C^{D) with respect to the norm 



1/2 



My= I \\u\\l2 + J2\\^Mh 



rl,2/ 



where the derivation YjU is intended in sense of distributions. For u G Wq ' {D; Y), we define 

the gradient Yu = {Yiu, . . . ,Ynu) and set \Yu\ = {\Yiu\^ -\ h |y„up)^^^ G L^iD). 

Setting Du := \Yu\ and H := W^' {D; Y), we define, for any VL <Z D, the energy Ef{Q) 
and the A:*'* eigenvalue Afe(O) of the operator Y^ + ■■ ■ + Y^, as in (|2.2|) and (|2.1[). The 



following existence result is a consequence of Theorem 2.4 



Theorem 2.12. Consider a bounded open set Z) C M*^ and a family Y = (Yi, . . . ,Yn) of 
C°° vector fields defined on an open neighbourhood U of the closure D of D. IfYi,...,Yn 
satisfy the Hormander condition on U, then for any A;GN, 0<c<|L'| and f G Lp'{D), the 
following shape optimization problems admit a solution: 



mm 



mm 



< Afc(i7) : i7 C -D, il quasi-open, \Q.\ < c\, 

< Ef{Q,) : n C D, Q. quasi-open, \Q\ < c}. 



(2.13) 
(2.14) 



1,2/ 



Proof. It is straightforward to check that the space H := Wq' {D;Y) and the apphcation 



Du := \Yu\ satisfy the assumptions of Theorem 2.4 The only non-trivial claim is the com- 
pact inclusion H C L'^{D), which follows since Yi,. . . ,Yn satisfy the Hormander condition 
on U. In fact, by the Hormander Theorem (see [17j). there is some e > and some constant 
C > such that for any <f G C^{D) 



where we set 



\v^\\h^ < C 



miH^ 



i</'IIl2 + ^\\Yj'^\\l2 

i=i 



(2.15) 



mo\\i + \^\'rdc 



1/2 



being (p the Fourier transform of (p. Let Hq[D) be the closure of C^{D) with respect to the 
norm || • \\h^- Since the inclusion L?'{D) C Hq{D) is compact, we have the conclusion. □ 



Remark 2.13. In the hypotheses of Theorem |2.12[ the following optimization problems have 
a solution: 

min < Afe(il) + A|r2| : Q C D, Q, quasi-open}, (2.16) 

min < Ef{Q) + A\Q\ : Q, C D, Q quasi-open}, 
where A; G N, A > and / G L^{D) are given. 



(2.17) 



Example 2.14. Consider a bounded open set D C M? and the vector fields X = -^ and 



Y 



' dy 



Since [X,Y] 



4-, we can apply Theorem 



2.12 



d_ 

dx 



and so, the shape optimization 



problem (2.17) has a solution Q C D. Assuming that iOs regular enough we may repeat 



the argument from Section 2.3 Indeed, suppose that y is a vector field on dfl and note 
that the map ^t = Id + tV is a differomorphism for t small enough. Defining Q^ = $t (Jl) 
and w the (strong) solution of 

-{dl + x^dl)w + 

where / G L?'{D), we have that 

dEfi^t) 



dt 



w = f, weW^^\n;X,Y), 
fw' dx, 



(2.18) 



t=o 



1 



(2.19) 



where w' is the weak solution of 

- {dl + x'^d'^) w' + w' = 0, w' + V-Vw£ W^'^{Qt; X, Y). 



Using (2.18) and integrating by parts in (2.19), we obtain 
dEfi^t) 



dt 



t=o 



Since 



1 

2 Jan 

d\nt\ 



dt 



{V ■ Vw) (n • {d^w,x'^dyw)) dn^. 

v-ndnK 



(2.20) 



(2.21) 



*=o Jan 

we have that the energy function w is a solution of the following overdetermined boundary 
value problem on the optimal set 



w 



[d1 + x^dy) w + w = f in il, 
= on dQ, 

2A on dn. 



(2.22) 



[ {n-{d.w,x^dyw))^^ 



an 



The characterization of the solutions of (2.22) is an open problem even in the case / = 1 



3 Spectral optimization for metric graphs 

In this section we study the problem of the optimization of the torsion rigidity of a one 
dimensional structure in M connecting a prescribed set of fixed points. Before we introduce 
the optimization problem we will examine some of the basic tools from the analysis of one 
dimensional sets. 

Consider a closed connected set C C M'^ of finite length 'H}{C) < cx), where by l-L^ we 
denote the one-dimensional Hausdorff measure in M . The natural choice of a distance on 
Cis 

dc{x,y)=\niU m)\dt : 7 : [0, 1] ^ M'^ Lipschitz, 7([0, 1]) C C, 7(0) = ^, 7(1) = ?/ 

which, in turn, gives a pointwise definition of a gradient 

I 71 / N ,. \u{v) — uix)\ 

\u'\{x) =limsupJ— ^^^^ ^, 

y^x a{x,y) 

which is a function in L'^{'H^), at least in the case when u : C — )• M is Lipschitz with respect 
to the distance dc- For any function u : C ^ M., Lipschitz with respect to the distance dc, 
we define the norm 

\hi(c) = I u^dV} + I \u'\^dV}, 

and the Sobolev space H^{C), as the closure of the Lipschitz functions on C with respect to 
this norm. By the Second Rectifiability Theorem (see [21 Theorem 4.4.8]) the set C consists 
of a countable family of injective arc-length parametrized Lipschitz curves 7^ : [0, k] — )• C, 
i G N, i.e. there is an ?^^ -negligible set iV C C such that C = NU (Uj 7i([0, /»])). On each 

curve 7j we have the chain rule -^u{^i{t)) = \u'\{'^i{t)) (see |8i Lemma 3.1] for a proof) 

and thus, we obtain the following expression for the norm of u S H^{C): 



\u\\ 



L"''"''^^!! 



rh 

|2 ~ " " 

' ' IC ~ JO 



|n(7.(t)) 



2 
dt. (3.1) 



Given a set of distinct points Di, . . . , D^ G M we define the admissible class Ac{Di, . . . , D^) 
as the family of closed connected sets C C M containing Di,...,Dfc. For any C G 
Ac{Di, . . . , Dfc) we consider the space of Sobolev functions which satisfy a Dirichlet condi- 
tion at the points Di: 

Hl{C; Di,...,Dk) = {ue H\C) : u{D,) = 0, j = 1 . . . , fc}. 

For the points Dj we use the term Dirichlet points. The Dirichlet Energy of the set C with 
respect to Di , . . . , D^ is defined as 

E{C]Di,...,Dk)= min - ( \u\^ dH^ - ( udU^ . (3.2) 

u&H^{C;Di,...,Dk) 2 Jc Jc 

We study the following shape optimization problem: 

min {E{C; Di,...,Dk) : C£ AciDi, ..., Dk), n\C) < l} . (3.3) 

Remark 3.1. We note that the admissible sets C can be reduced to the set of graphs 
embedded in M*^. For sake of simplicity, we limit ourselves to the case of three points 
Di,D2,D3 e R'^ (for the general result see [8j). Let C G Ac{Di,D2,D3) be such that 
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T-O'iC) < I and let r? : [0, a] — )• C be a geodesic in C connecting Di to D2 which we suppose 
that do not pass through D3. Let ^ : [0, 6] — )• C be a geodesic in C connecting D^ to Di 
and let ^3 G [0,6] be the smallest real number such that ^{h) G r/([0,a]). We define 

71 = ^|[o,«i], 72 = v{dc{Di,D2) - ■)\[o,i2], 73 = C|[o,y' 
where li and I2 are such that r/(/i) = ^(/s) and I2 = dciDi, D2) — h- 




Figure 1: The set C (on the left) and C (on the right). 

The curves 71 , 72 and 73 are geodesies in C which does not intersect each other in internal 
points (note that it is possible that one of them is degenerate, i.e. constant). Consider the 
set C = Ui7j([0,/i]) C C. By construction C' is connected and contains Di,D2 and D^. 
Let w e H^{C;Di,D2,D3) be a positive function and let v : [0,n^{C \ C')] -^ M be a 
monotone increasing function such that \{v < t}\ = ^{^{{w < r} fl L). By the Polya-Szego 
inequality (see [H Remark 2.6] or [H]), we have 



2 70 



hH^) 



\v'\ dx 



V dx < - 

2 ./r 



w'l'^dn^ 



wdU^. 



(3.4) 



Let a : \^,'H^{C \ C')] — )• W^ be an injective arc-length parametrized curve such that 
Im{a) n C" = cr(0) = x' , where x' £ C is the point where w\c' achieves its maximum. Then 
the closed connected set C = C"Uo"([0,'?^^(C\C")]) is admissible and has lower energy than 



C. In particular, in problem (3.3) with three fixed points, we can restrict our attention to 



sets, which are representations of metric graphs (i.e. combinatorial graphs with weighted 
edges) in W^. More precisely, we can consider graphs C such that 

1. C is a tree, i.e. it does not contain any closed loop; 

2. C has at most 6 vertices; if a vertex has degree three or more, we call it Kirchhoff 
point; 

3. there is at most one vertex of degree one for C which is not a Dirichlet point. In this 
vertex the energy function w satisfies Neumann boundary condition w' = Q and so we 
call it Neumann point. 

In the setting described above, the topology on the set of admissible graphs is quite 
natural, i.e. we say that C„ converges to C, if the weighted connectivity matrices of the 
graphs Cn converge to that of C, where the element rriij of the connectivity matrix M = 
is equal to the length of the edge connecting the two vertices Vi and Vj with the 



m,;,- 



'i])i] 



convention that rriij = +00 if the there is no edge connecting the two vertices and nnj = 0, 
if the two vertices coincide. It is quite clear that with this topology the set of connected 
metric trees of at most N vertices is compact. On the other hand, as the following example 
shows, the energy E{C,'D) is not semi-continuous. 
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Example 3.2. Consider the points Di = (0,0), D2 = (1,0) and D3 = (2,0) and the set 
Cn C M? consisting of the graphs of the functions y{x) = x{x — 1) for x G [0, 1] and 
Unix) = —^x{x — 2) for X G [0,2]. Passing to the hmit as n — ;■ 00, we have that the arc 
connecting Di to D3 passes through the Dirichlet point D2 which causes the energy to 
suddenly increase. 

Remark 3.3. The lack of semi-continuity does not necessarily imply the non-existence of 



a solution of (3.3), but suggests the nature of a possible counter-example. Following this 
idea, in [8j, was proved that if P = {Di, 02,0^} C M'^ is a set of points, with coordinates 
respectively (—1,0), (1,0) and (n, 0), and / = n -|- 2 is a given length, then, for n large 



enough, the problem (3.3) does not have a solution. 



In order to obtain an existence result for the problem (3.3), we consider, as in [8], 
in a larger class of admissible sets. Indeed, let F be a combinatorial graph with vertices 
{Vi}i=i^,,,^N and edges {eij}ij. We call F a metric graph, if to each edge Cjj is associated a 
positive real number lij which we interpret as the length of the edge. Thus, the total length 
of F is given by /(F) := Yli<j kj- 

A function u : F — )• M" on the metric graph F is a collection of functions Uij : [0, kj] — )• M, 
for 1 < i ^ j < N, such that: 

1. Uji{x) = Uij{lij — x), for each I < i ^ j < N, 

2. Uij{0) = UikiO), for all {i,j, A;} C {1, . . . , N}. 

We say that u is continuous (u G C'(F)), square integrable u G -^^(F) or Sobolev u G H^{T), 
if Uij is respectively continuous, square integrable or Sobolev on each edge eij. We also note 
that, if li G H^(T), then |u'| G i^(F) and so, we can define 

E{T;{Vi,...,Vk})= min \ I' \v! \'^ dV} - I' u dV} , (3.5) 

where Hq{T;{Vi, . . . ,Vk}) indicates the subspace of -ff^(F) of the functions vanishing on 
each of the vertices Vi,...,Vk and we also used the notation 

/ \u'\ d% := / / Wij\ dx, / udl-L ■=/ / Uijdx. 

Jv ^Jo Jv ^Jo 

We say that the continuous function 7 = {'~iij)i<ijLj<N : F — )• M'' is an im,mersion of the 
metric graph F into M , if for each 1 < i ^ j < N the function jij : [0, kj] — t- M is an 
injective arc-length parametrized curve. Given a set of distinct points Di, . . . , D^ G W^, we 
define the admissible set A{Di, . . . , Dk) as the set of metric graphs F for which there is an 
immersion 7 : F — )• M*^ such that 7(T^) = Di, where Vi, . . . ,Vk are vertices of F. In [8] the 
following result was proved. 

Theorem 3.4. Consider a set of distinct points Di, . . . , D^ G M"^ and a real number I such 
that there is a closed set C C M'^ which contains Di, . . . , Dk and such that TH}{C) < I. Then 
the following problem has a solution: 



mm 



{E{T;{Vi,...,Vk}) : F e A{Di, . . . ,Dk), /(F) < /}. (3.6) 



In some situations, we can use Theorem 3.4 to obtain an existence result for (3.3). 



Proposition 3.5. Suppose that Di, D2 and D^ be three distinct, non co-linear points in 
W^ and let I > be a real number such that there exists a closed set of length I connecting 



Di, D2 and D^. Then the problem (3.3) has a solution. 
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Proof. Let the graph F be a solution of (3.6) and let 7 : F — )• M be an immersion of F such 
that 'y{Vj) = Dj for j = 1, 2, 3. Note that if the immersion 7 is such that the set 7(F) C M'^ 



is represented by the same graph F, then 7(F) is a solution of (3.3) since we have 



E{r-{VuV2,V3}) = EiC;DuD2,Ds) 



Reasoning as in Remark 3.1, we can suppose that F is obtained by a tree F' with vertices 
Vi, V2 and V3 by attaching a new edge (with a new vertex in one of the extrema) to some 
vertex or edge of F'. Since we are free to choose the immersion of the new edge, we only 
need to show that we can choose 7 in order to have that the set 7(F') is represented by F'. 
On the other hand we have only two possibilities for F' and both of them can be seen as 
embedded graphs in M with vertices Di , D2 and D3 . □ 

Remark 3.6. Similarly to the existence proof of a classical optimal graph of Proposition 
above we believe that a more general result should hold: if Di, . . . , D^ are k distinct points 
in M'^ such that none of them can be expressed as a convex combination of the others, then 



(3.3) has a solution. We do not yet have a complete proof of this fact. 

"^ and let / > \Di — D2\ be a real 



Example 3.7. Let Di and D2 be two distinct points in 



number. Then the optimization problem (3.6) has a solution F which is actually a classical 



graph C given by the connected set (see Figure [2]) 



C = [Di,D2]U 



D1+D2 



,D3 



with 



D3 



D1+D2 



I- \Da 



D2\ 



V^ 



l-e 



9V4 



S 



— o- 
^3 



l-e 



Vo 



Figure 2: The optimal graph with two Dirichlet points. 

Example 3.8. Let Di, D2 and D^ be the vertices of an equilateral triangle of side 1 in 
i.e. 

/ \/3 \ _ /\/3 l\ ^ /V3 l\ 



D, 



,0 , i?2 



Z?.^ 



\ 3'/ \6'2/' " V^^/ 

We study in [8j the problem (lO) with V = {Di, L»2, -D3} and Z > \/3. We show that the 



solutions may have different qualitative properties for different / and that there is always a 
symmetry breaking phenomenon, i.e. the solutions do not have the same symmetries as the 
initial configuration D. Indeed, an explicit estimate of the energy shows that (see Figure [3]): 



1. if v3 < I < 1 + v3/2, we have that the solution of the problem (3.3) with D 
{Di, D2, D3} is of the form Fi; 



2. if / > 1 + \/3/2, then the solution of the problem (3.6) with V = {01,02,0^} is of 
the form F3. 
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Vi 





V2 Vi 

Figure 3: The optimal graphs for / < l+\/3/2, 1 = l+\/3/2, / > l+\/3/2 and I » l+\/3/2. 

4 Spectral optimization for Schrodinger operators 

Consider a bounded open set J7 C M"^ and a function / G L^($7). The Dirichlet energy 
related to a potential y > on Q is defined as 



Ef(V)= min - / \Vu\^ dx + - I u^V dx - / uf dx. 



[ v?Vdx- [ 
Jn Jn 



(4.1) 



A natural question, analogous to the problems considered in Section [2| is the optimization 
of Ef(V) under some integral constraint on V, of the form J V^ dx < 1. It is clear, from 
the definition of Ef, that for p > the minimum is achieved by K = 0. On the contrary, 
maximizing the energy under the same constraints gives the following results. 

• If p < la maximizing potential does not exist. In fact, for any p < 1, one may 
construct a sequence of functionals Vn such that J Vn dx = 1 and EfiVn) — )• 0, as 
n — )• 00. 



If p > 1 the optimal potential Vp exists and is given by 



K 



\u 



2/(p-l) . 



pP/(P-i)dx 



-1/p 



where u is the solution of the minimum problem 



mm - / |Vn| dx -\ — 
ueH^(n) 2 J^ 2 



u 



2P/{P-I)dx 



(p-l)/p 



uf dx, 



which is also the strong solution of —An + uVp = / in fi. 
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li p = 1 the optimal potential Vi exists and is given by 

where M = \\ui\\lo^(^q), u+ = {ui = M}, uj- = {ui = -M}, and ui e H^{0,)nH'^{n) 
is the unique minimizer of the functional Ji : L'^{Q) — t- M, defined as 

M'^) '■= n / \^u\'^dx+ -\\u\\1oom)- / ufdx. 
'^ Jn ^ Jn 

In particular, we have 

f dx — f dx = M, / > on lo^ , / < on uj^ . 

LJ^ J UJ — 

Example 4.1. Let Jl = (—1, 1) and / be a positive constant on Q. Then ui is positive 
and, by a symmetrization argument, it is also radially symmetric and decreasing. Thus, 
w+ = {—a, a) for some a G (0,1) and since |a;+|M = 1, we have that a = -^. Since 
n'(2jy) = and v!' = —f on (2^^^)^ ^^ have that (1 — 217)^/ — '^M, which uniquely 
determines M and so, the optimal potential Vi = -rr lf_^L j_\. 

J" I 2M ' 2M I 

When p < the minimization problem 

mmlEfiV) : y:O^[0,+oo], f V^ dx = l\ , (4.2) 

becomes meaningful. 

Proposition 4.2. Let i7 C M'^ &e a bounded open set and let f E L'^{Q). Then, for every 



p < 0, the problem (4.2) has a solution. 



Proof. By the definition of Ef(V), interchanging the two min operators, we find that the 
optimal potential Vp is given by 



Vp = |^|2/(P-i) . ( f \u\'P/(P~^) dx)''^' (4.3) 



-i/p 
In 
where u is the solution of the minimum problem 

If 1 / r \ (p~^)/p f 

min - / \Vu\'^dx + -{ / |u|2p/(P-i)d3; - / ufdx. (4.4) 

ueH^(n) 2 J^ 2 VJq / Jq 

Note that, since p < 0, the quantity q = 2p/{p — 1) is such that < g < 2. The existence 



of a solution for problem (4.4) is straightforward, which gives the existence of the optimal 
potential Vp through equality (4.3). D 



When we consider more general cost functionals FiV)^ like for instance spectral costs 
depending on the eigenvalues of the Schrodinger operator —A + V , the proof above cannot 
be repeated; nevertheless, using finer tools like 7-convergence for Dirichlet problems, the 
following more general result can be obtained (see [7]). 

Theorem 4.3. Consider a cost functional F : Bj^{^) — >■ M, where ;B+(J7) denotes the space 
of Borel measurable positive functions on il. Suppose that F is 

1. increasing, i.e. F[V) > F{W), whenever V > W; 
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2. lower semi- continuous with respect to strong convergence of the resolvents 

Rv = (-A + V)-^ : L\n) -^ L^{n). 

Then, for any p < 0, the optimization problem 

mmlFiV) : F:O^[0,+oo], f V^ dx = l\ , 

has a solution. 
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